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INTRODUCnON 


In  March  of  1964  Louisiana  State  University  Initiated  a  Nondestructive 
Testing  Institute.  This  Institute  was  formed  to  foster  Interest  In  academic 
programs  and  In  research  In  nondestructive  testing.  During  the  years  of 
1964  and  1965  the  Institute  was  mainly  concerned  with  conducting  short  courses 
for  civilian  and  service  personnel  In  the  field  of  nondestructive  testing. 

In  1966  application  was  made  to  the  THEMIS  program  for  a  grant  centering 
around  nondestructive  testing  technology,  One  of  the  projects  suggested  was 
entitled,  "Data  Analysis  and  Correlation  with  Digital  Computers".  A  portion 
of  the  project  statement  stated  that:  "Data  correlation  and  analysis  con¬ 
stitute  one  of  the  main  links  In  the  nondestructive  experiment.  Projects 
Involving  the  response  comparison  of  samples  to  a  standard  specimen  In  a 
vibrational  type  of  test,  or  In  a  test  by  means  of  the  reflection  of  ultra¬ 
sonic  waves,  or  In  the  back  scattering  test  of  nuclear  radiation  cannot  re¬ 
veal  meaningful  results  unless  proper  Interpretation  of  the  signals  from  the 
sensors  are  made.  These  analyses  can  be  efficiently  handled  by  high  speed 
digital  computers  using  methods  such  as  autocorrelation,  Fourier  transforms, 
and  such  statistical  tests  as  the  Chi  Square". 

The  alms  of  that  proposed  project  were  compatable  enough  with  the  alms 
of  the  ARFA  goals  so  as  to  be  of  Interest  to  that  agency.  Thus  the  mutua¬ 
lity  of  Interest  resulted  In  the  project  described  here. 

The  long  range  objective  of  the  project  was  to  develop  a  system  where¬ 
by  a  digital  computer  can  discern  the  difference  between  satisfactory  and 
unsatisfactory  basic  components  and  machine  parts. 

The  project  was  divided  Into  three  phases: 

Phase  I--Phase  I  provided  for  the  creation  of  a  computer  program  to 
analyze  and  display  data  when  steel  rods  were  subjected  to  various  energy 
exciting  media.  The  Investigation  apparatus  was  constructed  and  proven. 

Eight  (8)  months,  effective  September  12,  1968. 

Phase  II — The  computer  analysis  proved  to  be  dependable  In  reproduc¬ 
ibility  of  the  results  from  the  experimental  energy  envelopes.  The  computer 
program  was  flrmed-up  although  Improvements  In  the  program  continued  to  be 
made  throughout  the  life  of  the  contract. 

The  ability  of  the  system  to  differentiate  between  heat-treated  and 
nonheat-treated  samples  was  Investigated.  The  sensitivity  of  the  computer- 
program  to  determine  changes  In  dimensions,  structural  Integrity,  and  de¬ 
gree  of  stress  was  Initiated. 

Twelve  (12)  months,  effective  June  30,  1969. 

Phase  III--Investlgatlons  centered  primarily  around  the  determinations 
of  the  sensitivity  of  the  system  to  detect  dimensional  changes  and  degree 
of  stress  In  materials.  Structural  Integrity  was  Investigated.  The  results 
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of  hardness  tests  in  Phase  II  were  encouraging  enough  so  that  they  were  con¬ 
tinued  in  Phase  III.  Conversion  to  simulated  on-line  real  time  analysis 
was  completed  through  the  use  of  an  existing  hybrid  computer  system. 

Phase  III  extended  contract  for  additional  fourteen  (14)  months. 

In  addition  to  the  work  plan  described  above  the  researchers  inves+'i- 
gated  the  applicability  of  the  system  to  nondestructively  test  relatively 
long  (4  feet)  sections  of  thin  wall,  small  diameter  stainless  steel  tubing. 


Summary  of  Effect  and  Findings — Phase  I 

At  the  beginning  of  the  first  phase  of  the  effort  it  was  necessary  to 
design  and  construct  the  apparatus  and  instrumentation  to  accumulate  and 
store  data. 

Prior  to  building  the  apparatus  it  was  necessary  to  choose  the  types 
of  machine  elements  which  would  be  vibrated.  It  was  decided  to  use  Starrett 
tool  steel  machined  bars.  The  bars  were  all  18  inches  long  and  the  cross- 
sectional  dimensions  range  from  %  inch  by  ^  inch  to  ^  inch  square.  A  vise 
was  procured  to  hold  the  bars  while  they  were  being  vibrated.  In  order  to 
lessen  the  damping  characteristics  of  the  vise,  styrofoam  inserts  were  used 
between  the  vise  and  the  bars. 

The  next  problem  to  be  encountered  was  that  of  "ringing"  the  bars  at 
the  same  spot  with  the  same  magnitude  from  the  driving  device.  It  was  de¬ 
cided  to  use  a  solenoid  to  drive  a  threaded  pin  so  as  to  ring  the  specimens 
uniformly. 

To  record  the  signals  a  Computer  of  Average  Transients,  (CAT),  Model 
1000  manufactured  by  the  Technical  Measurements  Corporation,  was  obtained. 
The  CAT  Mode'  1000  has  a  time  corresponding  multi  channel  digital  memory 
whose  princl  >al  function  is  to  recover  repetitive  transient  signals  from 
random  background  noise. 

It  has  a  cathode  ray  tube  display  which  gives  a  visual  indication  of 
the  wave  pattern  at  the  end  of  the  accumulated  sweeps.  Additional  output 
of  paper  tape  printer  was  also  used. 

The  Technical  Measurement  Corporation  Model  500  Decimal  Printer  is  a 
combination  of  the  Hewlett  Packard  H51-565A  Printer  and  the  TMC  Control 
circuitry.  Basically,  the  Model  500  printer  prints  out  the  BCD  (binary 
coded  decimal)  coded  information  stored  in  the  CAT  unit.  The  printer  re¬ 
cords  the  complete  contents  of  the  CAT  memory,  when  directed  by  the  analyzer 
print  control.  This  permanent  record  is  printed  on  3-inch  standard  paper 
tape  at  the  rate  of  5  lines  per  second,  the  total  400  channel  memory  taking 
approximately  1  minute  and  20  seconds  to  print  out. 
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In  Figure  1,  A  indicates  the  DC  power  supply,  B  shows  the  bar, 

C  shows  the  Technical  Measurement  Corporation  CAT  unit,  D  Is  the  Gen¬ 
eral  Radio  Corporation  sound  level  meter,  E  Is  the  voltage  supply  for 
the  microphone  when  It  li?  being  used  as  a  pickup,  and  F  Is  the  Techni¬ 
cal  Measurement  Corporation  Model  500  Digital  Printer. 

In  summary,  the  work  accomplished  during  Phase  I  Includes  the 
checking  out  of  the  adaptability  of  the  Fast  Fourier  Transform  as  an 
adequate  analysis  tool  for  the  problem  under  consideration.  Various 
forms  of  frequency  energy  envelopes  were  Investigated  and  the  power 
spectral  representation  appears  to  offer  the  most  promising  signature 
for  comparison  purposes. 

The  computer  program  for  data  analysis  was  written  around  the  Fast 
Fourier  Transform  and  the  program  proved  effective  for  data  analysis. 

In  order  to  accomplish  reproducibility,  changes  had  to  be  made  In 
the  area  of  basic  experimental  set-up,  driving  medium,  recording  tech¬ 
nique,  and  record  length,  and  analysis  approach  In  establishing  a  firm 
signature. 


Summary  of  Effort  and  Flndlngs--Pha8e  II 

At  the  time  of  the  Initiation  of  Phase  II  a  primary  concern  was  to 
what  extent  the  pickup  system  and  the  analysis  program  could  discern  geo¬ 
metric  differences  In  tests  specimens.  However,  the  dimensions!  experi¬ 
ment  was  not  carried  out  Immediately  due  to  the  Interest  In  the  hardening 
determination  by  heat  treatment.  Although  the  goals  of  the  two  experi¬ 
ments  might  be  different  the  same  NDT  technique  was  applicable  since  the 
basic  principle  of  the  Investigation  was  based  on  the  frequency  deter¬ 
mination  of  a  vibrating  body.  The  body  vibrates  as  a  whole  In  one  of 
Its  natural  modes.  Many  physical  characteristics  of  the  body  may  be  deter¬ 
mined  from  the  characteristics  of  the  Induced  vibrations.  In  most  practi¬ 
cal  applications,  the  natural  frequencies  lie  largely  within  the  audible 
range,  from  a  few  to  about  20,000  cycles  per  second. 

When  a  thin  bar  Is  subjected  to  a  light  tapping  Impulse,  produced 
by  an  electronically  controlled  plunger.  It  response  Is  equivalent  to 
the  response  of  a  first  order  system  and  Is  Identical  with  Its  natural 
motion  from  an  Initial  position.  A  first-order  system  Is  one  whose 
motion  can  be  expressed  by  a  first-order  constant  coefficient  differen¬ 
tial  equation.  Since  the  Impulse  generates  a  very  short-time  Initial  con¬ 
dition,  the  response  produced  Is  the  natural  vibration  of  the  bar.  The 
resulting  natural  frequency  Is  sometimes  referred  to  as  the  resonant  fre¬ 
quency,  which  is  the  frequency  when  the  body  is  subjected  to  a  forced  vi¬ 
bration  whlth  the  same  frequency  as  Its  natural  frequency. 
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The  natural  frequency  of  mechanical  vibration  of  a  homogeneous  body 
is  controlled  by  a  nuinber  of  factors  which  may  be  generalized  by  2  factors 
in  the  expression: 

Frequency  =  (shape  factor)  x  (physical-constants  factor) 

Shape  factor  Includes  the  geometric  dfitiign  of  the  body  and  the  dimen¬ 
sional  factors  of  length  and  thickness.  The  physical-constants  factor  in¬ 
cludes  the  modulus  of  elasticity,  density  and  Poisson's  ratio  for  the  test 
material. 

It  is  axiomatic  that  two  bodies,  identical  in  every  respect,  will  mani¬ 
fest  identical  natural  frequency  in  vibration.  If  the  shape  factor  of  a 
number  of  bodies  is  held  constant,  the  natural  frequency  of  vibration  of 
each  of  the  bodies  will  be  a  discrete  n^asure  of  the  physical-constants 
factor.  On  the  other  hand,  if  the  effect  of  the  shape  factor  on  frequency 
can  be  determined,  compensation  for  differencies  in  shape  or  size  may  be 
applied. 

Every  object  has  many  modes  of  vibration.  In  addition  to  the  funda¬ 
mental  mode  it  maybe  expected  that  the  complexity  of  modes  and  overtones 
would  coii4>licate  the  resonance  vibration  test.  However,  it  is  not  the 
case  in  this  experiment  because  the  impulse  excitation  would  not  excite  the 
higher  order  overtones  at  a  simply  supported  configuration.  The  funda¬ 
mental  frequency  of  vibration  is  usually  the  most  conveniently  detectable 
frequency. 

Vibration  testing  was  used  in  flaw  detection  by  others.  McMasters^ 
reported  the  example  of  the  flaw  simulated  by  a  lateral  saw  cut  in  the 
center  of  a  square  bar.  The  method  of  detecting  the  flaw  consisted  of 
measuring  the  response  frequencies  of  the  fundamental  mode  after  placing 
the  bar  into  separate  orientations.  Table  1  was  prepared  from  actual 
measurements  of  a  simulated  and  progressively  deepened  flaw. 

It  can  be  seen  that  there  is  a  difference  in  the  response  frequency  as 
the  depth  of  the  flaw  Increases. 

The  hardening  of  a  steel  specimen  by  heat  treatment  could  be  thought  of 
as  a  body  experiencing  uniform  microscopic  flaws.  When  steel  is  in  the  non- 
hardened  conditions  it  is  essentially  pearllte.  This  microstructure  is  a 
laminar  mixture  of  ferrite  and  carbide.  The  lighter  color  is  the  ferrite 
and  the  darker  color  is  carbide.  There  is  nearly  127,  carbide  and  slightly 
more  than  88%  ferrite  in  the  mixture.  Since  the  carbide  and  ferrite  form 
simultaneously,  they  are  intimately  mixed.  Characteristically  the  mixture 
is  laminar.  That  is,  it  is  composed  of  alternate  layers  of  ferrite  and  car¬ 
bide.  The  resulting  microstructure  is  called  austenite. 


^McMaster,  Robert  C. ,  Nondestructive  Testing  Handbook,  Vol.  2, 
page  51-12. 
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'  TABLE  1 

I 

'Frequency  In  Hertz  with  Orientations*  >  ,  , 


Depth  of 

Flaw,,  in.  , 

A 

I 

B 

Frequency 

1  Difference 

0  '1  ? 

1552 

1552 

0  ' 

1/16 

1533 

'  1541 

'  8 

1/8 

1A69 

1533 

64 

3/16  ,  ,  .  ' 

1400 

1513 

113 

,1/A'  '  , 

1327 

'  1507  1 

180  , 

*Orientatipns  of  the  same  bar  are 
Position  B  is  vertical. 

different.  Position  A, Is 

horizontal. 

After  heat  treatment,  as  temperature  decreases  thie  free  energy  difference 
between  the  austenite  aftd  the  combination  of  ferrite  plus  carbide  becomes 
I  larger.  Consequently,  the  austenite  is  mofe  unstable.  At  sufficiently  low 
temperatures,  but  usually  still  above,  room  temperatures the  face-centered 
structure  does  change  to  a  body-centered  structure,  but  all  of  the  carbon  pre¬ 
sent  remains  in  splutlon.  Consequently,  thfe  resulting  body-centered  structure 
is  tetragonal  rdther  than  cubic.  It  is  then  called  martenSitej  i 

'  I 

1  •  'I 

Since  martensite  has  a  noncuble  structure  .it  is  reasonable  to  expect  that 

'  any  energy  signature  bbtained  from  a  cubic  structure  then  noncuble  structure 
would  have  a  differential.  It  is  believed  that  this  difference  in  structure 
results  in  the  difference  Ln  frequency.  Until  further  work  is  done  it  caii  ' 
not  be  argued  that  the  ffoqueiicy  it  entirely  attributable  to  grain  structure. 

It  might  be  attributable  to  the  resulting  stresses  set  up'  in  the  steel.  ' 

The  first  three  months  of  the  Phase  II  were  spent  in  continuing  develop¬ 
ment  of  the  system  procedure  to  differentiate  between  very  small  hhanges  in 
configuration  and  integrity  of  the  steel  bars  that  were  being  tested.  During 
that  time  period  the  ability  to  reproduce  results i  when  the  same  bar  is  vi - 
I  brated  after  being  removed  from  the  holder  and  reinserted  in  the  holder,  was 
perfected.  ,  '  ,  i 

I  'll 

Following  that,  investigation  into  the  ability' of  the  system  to  determine 
a  frequency  correlation  of  a  sample  after  being  heat  treated  was  initiated, 
fhe  first  step  consisted  df  choosing  two  bars  and  subjecting  them  to  the  com¬ 
puter  analysis  both  before  and  aftet  heat  treating.  The  results  offered 
enough  promise  to  pursue  the  problem.  The  second  .step  consisted  of  fabricating 
eight  sample  bars,  each  eighteen  inches  long  from  the  same  piece  of  one-half 
inch  bar  stock;  ,  ' 

The  bars  were  subjected  to  a  Brinell  hardness  test  to  determine  their 
before  heat  treated  hardness.  Energy  signatures  were  then  made  of  each  bar 
and  the  bars  were  then  heat  treated.  After  heat  treating,  the  hardness  of 
the  bars  was  again  determined  and  the  energy  signatures  again  r^ecorded.  The 
bars  were  then  Xfrayed  to  determine  if  there  were  any  'discontinuities  in- 
troduced  as  a  result  of  the  heat  treating  operation.  Table  2  shows  the  re¬ 
sults  of  these  determinations.  ' 

It.' 

Ill/  ’  :  I  ' 
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' TABLE  2 

Sample  Hardness  and  Frequency- -Run  One 


Rod 

Identi 

Brlnell  ' 
before  H.T. 

Brlnell 
after  H.T. 

Response 
freq,  before 
H.T.,  Hertz 

Response 
freq.  after 
H.T. , Hertz 

F 

'  57.7 

56.2  , 

588 

584 

G 

,  61.0 

59.0 

588 

580 

k 

58.4 

60.3 

592 

588 

1  I 

,  59.7  , 

61.0 

592 

580 

J  , 

57.8 

<  54.5 

592 

576 

L 

58.4 

58.8 

588 

580 

M 

59.2 

'  55.1 

592 

572 

N 

59.7 

58.0 

I 

592 

584 

It  will  be  realized  that  the  magnitude  of  the  hardness  differences  en¬ 
countered  In  the  sample  were  slight.  The  trend  of  the  tests  suggested  that 
some  correlation  between  hardness  and  frequency  response  of  the  samples. 
Little  credence  was  given  to  these  results  Alnce  certain  errors  became  appa¬ 
rent  after  the  runs  were  made.  i 

I 

Therefore,'  In  order  to  verify  the  relationship,  another  set  of  samples 
was  prepared.  This  time  the  samples  were  made  from  1095  stock,  a  steel 
more  conducive  to  heat  treating  than  the  1020  used  In  run  one.  These  are 
prefixed  by  the  letter  "A".  Figure  2  describes  the  data  collected  during 
this  experiment. 

{  I 

'  '  TABLE  3 

I  ; 

'  Sample  Hardness  and  Frequency — Run  Two 


Rod 

Ident. 

Response 
freq.  before 
H.T., Hertz 

Response 
freq.  after 
H.T.,  Hertz 

Rockwell  C, 
before  H.T. 

Rockwell  C, 
after  H.T. 

.  AA 

945.3  ' 

937.5 

8.6 

42.6 

AB 

960.9 

947.3 

8.6 

36.5 

AC 

'  957.0 

945.3 

8.6 

36.5 

AE 

953.1 

945.3  ' 

8.6 

27.3 

AG 

960.9 

945.3 

8.6 

47.5 

AH  1 

946.3  . 

921.9 

8.6 

54.0 

AI 

968.8 

953.1 

8.6 

33.1 

AJ 

945.3 

937.5 

8.6 

36.9 

AK 

945.3 

937.5 

8.6 

6.0 

AL 

937.5 

937.5 

8.6 

2.2 
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Table  4  reduces  the  data  In  Table  3  to  a  differential  hardness  related 
to  a  differential  frequency  and  this  Information  Is  plotted  on  Figure  2. 

Figure  2  Indicates  the  data  collected.  The  degree  of  relationship  be¬ 
tween  the  two  above  mentioned  variables  Is  shown  below. 

The  data  shown  In  Figure  2  was  obtained  by  taking  10  bars  and  subjecting 
these  bars  to  various  degrees  and  Intensities  of  heat  treatment. 

TABLE  4 

Differential  Hardness  vs.  Differential 
_ FrecjUency — Run  Two _ 


^  Hardness,  Rockwell  G,  A  Frequency,  Hertz 

Rod  rfter-before  after-before 


AA 

34.0 

-7.8 

AB 

27.9 

-3.6 

AC 

27.9 

-11.7 

AE 

18.7 

-8.8 

AG 

38.9 

-15.6 

AH 

45.4 

-24.4 

A1 

24.5 

-15.7 

AJ 

28.3 

-7.8 

AK 

-2.6 

-7.8 

AL 

-6.4 

-0.0 

The  data  were  subjected  to  a  conventional  correlation-regression  treat¬ 
ment  to  obtain  a  mathematical  model.  The  model  Is: 


y^  =  -4.11  -  .30x 


with  a  coefficient  of  correlation  of  .78,  where  the  variables  are  as  defined 
previously. 


The  analysis  of  variance  of  regression  Is  as  follows: 


MKlMR 


—  w 


TABLE  5 
Analysis  of  Variance 


Source 

Degrees 
of  Freedom 

Sum  of 
Squares 

Mean 

Square 

F 

Ratio 

Regression 

1 

232.3 

232.3 

12.1 

Residual 

8 

153.3 

19.2 

Total 

9 

385.6 

During  Phase  11  It  became  apparent  that,  for  the  data  analysis  system 
to  be  of  any  practical  use,  a  faster  system  should  be  sought.  Coincidental 
with  this  realization,  the  LSU  hybrid  computer  system  became  operational. 

This  computer,  supported  by  the  Department  of  Defense,  under  a  THEMIS  grant. 
Involved  an  expenditure  of  over  $750,000  for  hardware  alone. 

A  hybrid  system  receives  analog  Information  such  as  voltages  which  are. 

In  this  case,  proportional  to  the  sound  signature.  The  LSU  hybrid  system 
uses  a  'ilgma  5  for  the  digital  component  and  an  EAI  680  for  the  analog  com¬ 
ponent.  Using  the  hybrid  system.  It  Is  necessary  merely  to  record  on  magnetic 
tape,  one  cycle  of  an  energy  signature.  The  tape  Is  then  used  as  Input  to 
the  hybrid  system.  Figure  3b  depicts.  In  block  diagram  form,  the  schen;atlc 
of  the  hybrid  analysis.  The  digitizing  Interval  can  be  varied  by  a  simple 
change  of  circuit  elements.  The  play  back  speed  can  also  be  varied.  The 
analog  signal  Is  converted  to  digital  Information  through  an  A  to  D  converter 
(ADC).  This  Information  Is  In  turn  stored  sequentially  in  the  memory  core  of 
the  digital  part  of  the  hybrid  computer.  As  soon  as  a  preset  number  of 
sampling  points  are  obtained,  the  calculations  begin.  The  final  result  can 
be  either  display  on  an  oscilloscope,  or  on  a  x-y  plotter,  or  on  a  print 
plot.  The  print  plot  was  chosen  because  of  its  simplicity  and  clarity  of 
locating  the  resonance  peak.  It  Is  also  the  fastest  way  to  obtain  a  permanent 
record. 

In  an  attempt  to  verify  the  computer  program  for  the  hybrid  system,  the 
signature  of  a  pulser  generated  frequency  (625  Hertz)  was  run  on  the  hybrid 
system.  A  portion  of  the  results  are  shown  on  Figure  4.  It  has  been  shortened 
for  presentation  purposes  and  the  line  connecting  the  ’’P's"  has  been  drawn 
in  for  clarification.  This  printout  is  analogous  to  that  obtained  from  the 
calcomp  x-y  plotter  used  with  the  previous  system. 

In  Figure  4  the  N  is  the  number  of  digitized  sampling  points  of  the 
sound  signature.  LG  (LAG)  is  the  autocorrelation  lag.  TIME  RES.  is  the  time 
interval  between  sampling  points  in  milliseconds.  The  Nyquist  frequency  is 
the  maximum  frequency  that  can  be  observed  with  the  given  time  resolution. 

FREQ.  RES.  is  the  frequency  resolution  interval  between  consecutive  points 
in  the  power  spectrum.  The  R  column  is  the  list  of  normalized  autocorrela¬ 
tion  function  values.  The  P  column  is  the  list  of  normalized  power  spectral 
values. 
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FIGURE  3 

BLOCK  DIAGRAMS  OF  COMPUTER  AN^YSIS  CONFIGURATION 


I  Computer  of 


Average 

Transients 


Paper 
I  Tape 
I  Printer 


i 

i 

I  .  .  Key 
I  Punch 


(a) 


Original  Method 


(b)  . 

Hybrid  Configuration 


II 


FIGURE  4 


The  ratio  of  P/ (P  SCALE)  (see  heading  for  P  SCALE)  gives  the  actual  P 
value  calculated.  The  ratio  of  R/(R  SCALE)  (see  heading  for  R  SCALE) 
gives  the  actual  R  value  calculated.  The  *  indicates  P  and  R  coincide. 
The  autocorrelation  function  is  used  to  indicate  the  degree  of  constancy 
of  a  particular  signal  quantity  as  a  function  of  time.  The  power-density 
spectrum,  which  represents  a  measure  of  the  distribution  of  energy  in  the 
frequency  spectrum,  is  the  Fourier  transform  of  the  autocorrelation  func¬ 
tion. 


The  programming  and  theoretical  considerations  resulting  from  Phase  II 
are  included  in  Appendix  B. 

A  number  of  conclusions  were  drawn  from  the  work  performed  during  Phase 
II. 

1.  There  appeared  to  be  a  predictable  relationship  between  changes 
in  the  hardness  of  certain  carbon  steels  and  changes  in  response 
frequencies. 

2.  Hybrid  computational  equipment  can  give  on-line  analyses  of  these 
situations. 

3.  Frequency  changes  could  be  attributable  to  differences  in  hardness 
and  stress  differences. 

In  any  experimental  investigation  such  as  that  described  here,  a  mathe¬ 
matical  model  is  sought  to  define  the  cause  and  effect  aspects.  Ideally, 
if  a  coefficient  of  correlation  of  1.0  would  be  obtained,  it  could  be  assumed 
that  only  a  change  in  hardness  causes  a  change  in  response  frequency.  In 
run  two,  a  coefficient  of  correlation  of  .78  indicates  that,  to  some  small 
degree,  other  factors  are  present.  One  of  these  factors  might  have  been  the 
sample  size.  The  other  might  have  been  the  fact  that  a  variable  or  variables 
other  than  hardness  is  causing  a  change  in  frequency.  But,  in  view  of  the 
coefficient  of  correlation  obtained,  the  majority  of  the  variation  is  probably 
accounted  for  by  a  change  in  hardness  alone. 

The  procedure  used  for  analysis  during  Phase  II  consisted  of  collecting 
the  energy  signatures,  averaging  a  number  of  them,  and  then  digitizing  them 
on  paper  tape.  Cards  were  then  cut  from  the  paper  printout  tape  and  the  data 
run  on  the  program  using  the  IBM  360  Model  65  system.  This  was  a  time  con¬ 
suming  process  requiring  averaging  3  days  from  data  acquisition  to  receipt  of 
plots.  With  the  program  prepared  during  Phase  II,  every  indication  pointed 
to  the  successful  adaption  to  the  existing  LSU  hybrid  computer  to  analyze 
the  data.  It  was  still,  however,  a  simulated  on-line  system  since  a  sound 
recording  tape  was  used  to  record  the  sound  and  it  had  to  be  transferred  to 
the  hybrid  system  to  be  analyzed.  If  tape  could  be  eliminated,  the  direct 
sound  signature  results  could  be  obtained  immediately.  This  would  change 
from  a  simulated  on-line  system  to  a  truly  on-line  system,  that  is,  results 
could  be  obtained  within  one  minute  or  less  after  vibration  of  the  sample. 


13 


Following  the  frequency-hardness  investigations,  the  frequency-dimension 
project  was  started  during  Phase  II. 

Fourteen  bars  were  cut  and  end  center  drilled  for  a  lathe  mounting. 

The  diameter  of  each  bar  was  reduced  incremei tally  to  ascertain  the  sensitivity 
of  the  analysis  procedure  to  react  to  dimensional  differences.  These  runs 
were  made  on  the  hybrid  system. 

Fourteen  rectangular  bars  were  prepared  to  stress  using  a  tensile  testing 
machine.  Before  and  after  signatures  were  made  to  determine  whether  or  not 
the  procedure,  as  developed,  could  differentiate  between  the  before  and  after 
case  and  the  nature  of  that  relationship.  Results  of  this  work  are  described 
in  the  next  section. 


i 
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TECHNICAL  DISCUSSI(»}  OF  WORK  PERFORMED 


PHASE  III 


During  Phase  III,  the  following  topics  were  Investigated, 

1.  The  stress  analysis  was  completed  on  twenty  rectangular  bars 
(six  1010  carbon  steel  and  fourteen  1040  carbon  steel).  The 
results  Indicate  that  there  Is  a  promising  relationship  between 
the  stress  and  resulting  frequency  of  the  bars  and  that  this  re¬ 
lationship  Is  only  negligibly  dependent  on  the  dimensional 
changes  of  cross-sectional  area  and  the  bar  length  that  result 
from  the  stress. 

2.  A  non-stress  analysis  of  the  dimensional  property  of  length  was 
performed  on  a  carbon  steel  round  bar  to  ascertain  the  effect 
of  a  change  In  length  while  the  dlan^ter  remained  unchanged. 
During  this  period  a  frequency  resolution  Improvement  was  made 
to  the  computer  program  to  give  more  sensitive  results.  This 
fact  added  to  the  attractiveness  of  the  acoustic  technique. 

3.  An  analysis  of  the  dimensional  properties  of  the  diameter  of 

a  non-stress  carbon  steel  round  bar  was  performed  to  ascertain 
the  effect  of  a  change  In  diameter  while  the  length  remained 
unchanged . 

4.  Rolled  steel  specimens  were  tested  to  ascertain  If  the  hardness 
could  be  determined  by  the  acoustical  technique.  These  specimens 
were  received  from  Frankford  Arsenal. 

5.  An  attempt  to  use  the  acoustical  technique  to  monitor  the  Inte¬ 
grity  of  long,  thin  wall,  small  diameter  tubing  was  made.  The 
stainless  steel  tubes  were  sent  from  WADCO  (a  subsidiary  of 
Westlngh3use  Electric  Corporation)  In  Richland,  Washington. 

6.  A  1/60"  hole  was  drilled  through  the  diameter  center  axis  of 
two  round,  carbon  steel  bars  to  determine  If  a  frequency  change 
would  result  between  them  and  Integral  bars. 


Stress  Analysis 

The  stress  analysis  Investigation  was  Initiated  by  loading  a  1010 
carbon  steel  rectangular  bar  on  the  Instrom  tester  model  TT  type  D.  Fig- 
urge  5  shows  this  experimental  configuration.  The  bar  was  loaded  by 
applying  Increasing  tensile  stresses  In  percentages  of  Its  yielding  stress. 
At  each  Increment  the  bar  was  struck  by  the  solenoid  and  the  resulting  sound 
was  recorded  on  the  tape  recorder.  This  recording  was  then  Introduced  Into 
the  hybrid  computer  to  be  analyzed.  The  response  frequency  of  the  bar  was 
determined  and  printed.  Six  1010  carbon  steel,  rectangular  bars  were  tested 
In  two  groups  of  three.  The  results  of  group  one  are  listed  In  Table  6 
with  the  resulting  curve  In  Figure  6.  Once  the  configuration  for  group  one 
was  established,  group  two  was  then  run  with  more  attention  paid  to  detail. 
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Figure  5 

Configuration  Using  Instron  Testing  Machine 
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This  data  Is  shown  In  Table  8.  A  curve  of  the  results  Is  shown  in  Figure  7. 
Note  the  relation  between  the  experimental  and  theoretical  frequencies. 

The  theoretical  model  Is  developed  in  Appendix  C. 


TABLE  6 

LEAST  SQUARES  ANALYSIS  FOR  DATA  ON  BARS 
SHOWN  IN  TABLE  8 


y  =  7.  Yield 


X  »  Exper.  Freq. 

X  =  Theoretical  Freq, 

Bar  I 

y  -  143.5  +  .588  x 

R  =  .99 

y  -  126.18  +  1.59  x 

R  «  .99 

Bar  II 

y  “  133.6  +  .63  x  • 

R  «=  .99 

y  -  124.3  +  1.72  x 

R  -  .996 

Bar  III 

y  =  138.2  +  .49  X 

R  «  .98 

y  »  122.6  +  1.75  x 

R  -  .99 
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TABLE  8 


EXPERIMENTAL  VS  THEORETICAL  FREQUENCIES -STRESS 

Bar  I 
460  g 


Exper. 


Theoretical 


Yield 

Load 

Width 

Depth 

Length 

Freq. 

Freq. 

0 

.9897 

.2518 

11.324 

117 

II  II 

5 

500 

.9887 

.2518 

11.326 

II  II 

131 

11 

1020 

.9885 

.2515 

11.326 

148 

•  !-» 

V£) 

144 

VO  ro 

VO  ov 

23 

2050 

.9887 

.2518 

11.330 

156 

VO  U> 

167 

34 

3050 

.9890 

.2515 

11.334 

164 

186 

00 

44 

4000 

.9885 

.2516 

11.335 

171 

■T 

203 

+ 

57 

5100 

.9887 

.2517 

11.337 

179 

.588 

221 

I-* 

68 

6100 

.9890 

.2515 

11.340 

187 

235 

cn 

80 

7130 

.9888 

.2514 

11.346 

187 

250 

90 

8040 

.  9895 

.2515 

11.350 

195 

262 

Bar  II 

447  g 

0 

25 

1.002 

.2521 

11.682 

132 

<< 

119 

71  Vj 

10 

900 

1.002 

.2521 

11.687 

140 

II 

142 

II  II 

20 

1820 

1.0018 

.2520 

11.690 

148 

I-* 

CO 

162 

•  K’ 

VO  ho 

30 

2700 

1.0020 

.2520 

11.693 

156 

u> 

• 

181 

'O 

o>  • 

40 

3600 

1.0015 

.2519 

11.698 

156 

<J* 

197 

CO 

50 

4500 

1.0018 

.2519 

11.700 

164 

212 

+ 

60 

5400 

1.0015 

.2519 

11.700 

171 

• 

CV 

227 

• 

70 

6300 

1.0015 

.2519 

11.702 

179 

U) 

X 

240 

ro 

X 

Bar  III 

455  g 

0 

100 

.9982 

.2521 

11.796 

II  II 

119 

II  II 

10 

900 

.9980 

.2520 

11.800 

140 

• 

140 

•  !-» 

20 

1800 

.9980 

.2520 

11.808 

148 

VO  CO 

00  00 

160 

vO  Ni 
VO  ro 

30 

2700 

.9975 

.2520 

11.812 

156 

• 

ro 

178 

<j* 

40 

3600 

.9970 

.2520 

11.814 

156 

+ 

195 

+ 

50 

4500 

.9965 

.2519 

11.815 

164 

• 

210 

I-* 

60 

5400 

.9965 

.2519 

11.818 

171 

•P> 

VO 

224 

70 

6300 

.9965 

.2519 

11.822 

171 

X 

238 

80 

7200 

.9962 

.2519 

11.829 

179 

251 

X 

90 

8100 

.9962 

.2519 

11.836 

179 

263 

y  =  %  yield;  x  =  exper.  freq.;  x  =  theoret.  freq. 
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7.  .Yield 


Fourteen  1040  stainless,  carbon,  rectangular  steel  bars  were  then  run 
to  determine  whether  the  frequency  change  was  due  to  a  change  in  cross- 
sectional  area  or  change  in  stress.  The  bars  were  tested  in  the  same 
manner  as  the  1010  carbon  steel  bars  referenced  to  previously.  However, 
because  they  were  essentially  strainless,  their  cross-sectional  area 
had  little  (.3%)  measurable  change.  Due  to  this  property  of  the  bars, 
the  results  Indicate  that  the  relationship  is  mainly  between  stress  and 
response  frequency  of  the  bars.  The  data  and  results  for  these  bars  are 
listed  in  Table  9  with  an  accompanying  curve  in  Figure  8.  The  complete 
data  on  this  run  is  Included  in  Appendix  D. 
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RESPONSE  FREQUENCY  FOR  1040  BARS  UNDER  VARIOUS  LOADS 
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Dimensional  Study 


A.  Length  Variation 

In  search  for  the  correlation  between  the  length  of  a  bar  and  Its  response 
frequency,  a  round  carbon  steel  bar,  approximately  12"  long  was  prepared.  It 
diameter  was  constant.  The  bar  was  freely  supported  at  two  ends  with 
styrofoam  In  contact  at  each  end.  There  was  no  rigid  restrain  at  the  sup¬ 
porting  ends.  The  weight  of  the  bar  and  the  friction  at  the  ends  are  the  only 
means  to  keeping  the  bar  In  position  vdien  struck  from  the  side  at  Its  mid¬ 
point  by  the  solenoid  plunger.  The  energy  envelope  of  the  bar  was  recorded 
and  analyzed.  The  frequency  of  the  bar  at  each  decremental  length  was  ob¬ 
tained  and  the  results  tabulated  In  Table  10.  The  relationship  between  the 
experimental  frequency  and  the  length  appear  In  good  agreement  with  the  theore¬ 
tical  model.  The  theoretical  response  frequency  of  a  round  bar  with  hinged 
ends  Is  (see  )^pendlx  E). 


4 

In  which,  I  B  ttD  /64  for  round  bars.  The  length  and  frequency  relationship 
becomes 


f 


l6(L/n)^ 


(2) 


If  all  terms  are  held  constant  except  the  length  and  the  frequency,  their 
relationship  Is 


The  experiment  result  verified  this  Inverse  square  relationship.  However, 
the  magnitude  of  the  experimental  and  theoretical  frequencies  disagreed  by 
a  constant  factor.  Both  theoretical  and  experimental  frequencies  plotted 
as  a  function  of  length  are  shown  on  Figure  10.  Further  study  of  this  dis¬ 
crepancy  reveals  that  the  factor  n  In  Equation  (2)  played  an  Important 
role  In  the  analysis,  n  Indicates  whether  the  bar  vibrates  In  a  fundamental 
mode  with  n  =  1,  or  In  the  second  harmonic  with  n  2,  or  In  higher  order 
harmonics  with  n  >  2.  When  n  a  i,  the  vibration  has  nodes  at  0  and  L.  This 
Is  a  hlnged-end  model.  When  n  =  2,  the  nodes  are  at  0,  L/2  and  L.  The  ex¬ 
perimental  result  Indicated  that  the  n  value  was  between  1  and  2.  The  reason 
for  the  percullar  node  positions  resulted  from  the  bar  sopports.  They  were 
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not  hinged  ends.  The  bar  ends  were  free  to  vibrate  on  the  surface  of  the 
styrofoam  when  the  bar  was  struck  at  midpoint  from  the  side.  The  distance 
between  the  nodes  was  found  to  be  0.662  L.  This  means  n  «=  1.5186.  With 
this  new  value  of  n,  the  theoretical  curve  coincided  with  the  experimental 
curve.  These  are  illustrated  in  Figures  9  and  10.  For  this  reason,  if 
the  experimental  set  up  was  a  hinged  end  configuration,  the  experimental 
frequencies  would  agree  with  the  n  =  1  theoretical  values. 


TABLE  10 


CHANGE  IN  LENGTH  OF  ROUND  BAR 


DIAMETER  0.7000" 


X 

Length 

In 

Weight 

Grams 

^2 

Theoretical 

^1 

Frequency,  Hz 
Experimental 

12.043 

595 

384 

883 

11.945 

590 

390 

898 

11.844 

585 

397 

932 

11.736 

580 

404 

942 

11.626 

575 

412 

951 

11.535 

570 

418 

966 

11.438 

565 

425 

981 

11.315 

560 

434 

996 

11.186 

555 

444 

1020 

11.064 

545 

455 

1035 

Least  Squares  Analysis  with  a  linear  model: 
Experimental  y^  =  2700.3  -  150.34  x 
R  =  -.99 

Theoretical  =  1245.2  -  71.62  x 
R  =  -.997 


During  this  period  of  the  experiment  the  computer  program  sensitivity 
has  been  improved.  It  is  known  that  the  frequency  resolution  limitation  had 
an  effect  on  the  sensitivity  of  obtaining  the  experimental  frequencies.  The 
experimental  frequency  calculations  are  based  on  the  power  spectural  analysis 
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of  the  time-series  data.  The  final  frequency  resolution,  which  Is  the 
smallest  frequency  change  that  can  be  resolved,  Is  a  function  of  sampling 
Interval  of  the  analog  data  and  the  maximum  lag  of  the  autocorrelation 
function.  If  the  sampling  Interval  Is  made  large,  the  Nyqulst  frequency 
will  be  lowered.  However,  the  frequency  resolution  will  be  Improved.  If 
the  response  frequency  of  the  bar  continue  to  fall  within  the  lowered 
Nyqulst  frequency,  It  appears  wise  to  Increase  the  sampling  Interval  for 
the  analog-to-dlgltal  conversion.  The  other  approach  was  to  Increase 
the  maximum  lag  number  of  the  autocorrelation  function.  However,  since 
the  Hybrid  computer  has  a  limited  storage  location  (32  K  memory),  an  In¬ 
crease  of  the  maximum  lag  from  256  to  512  will  over  fill  the  core  storage. 
Since  a  256  to  512  Increase  Is  a  minimum  step  Increase,  the  alternative 
remedy  Is  to  reduce  the  sampling  record.  Three  thousand  points  of  Infor¬ 
mation  were  gathered  for  each  run.  Since  the  nature  of  the  data  Is  a 
transient  one,  1000  to  2000  points  may  suffice.  A  modification  of  the 
computer  program  has  achieved  this  goal  and  Is  reproduced  as  Appendix 
F.  It  was  found  that  by  Increasing  the  sampling  Intervals  from  0.4  msec 
to  0.5  msec,  a  reduction  of  frequency  resolution  from  8.9  Hz  to  7.8  Hz 
was  produced.  Further  reduction  was  made  by  revising  the  computer  pro¬ 
gram  to  give  a  greater  max  lag  In  correlation  analysis.  The  result  gave 
a  3.9  Hz  resolution.  However,  the  Nyqulst  frequency  was  limited  to  1000 
Hz.  The  summary  of  the  sampling  time,  max  lag  and  frequency  resolution 
Is  shown  as  follows: 


Sampling 

Time  (msec) 

Maximum 

Lag 

Frequency 
Resolution  (Hz) 

Nyqulst 
Frequency  (Hz) 

.4 

256 

9.7656 

1250 

.4 

512 

4.8828 

1250 

.5 

256 

7.8125 

1000 

.5 

512 

3.90625 

1000 

B.  Diameter  Variation 

Similar  efforts  have  been  made  for  the  correlation  between  the  diameter 
of  a  bar  and  Its  response  frequency.  A  carbon  steel  round  bar  was  prepared 
for  the  experiment  by  altering  its  diameter  on  a  lathe  while  the  length  re¬ 
mained  constant.  The  diameter  was  cut  in  decrements  from  0.7078"  to  0.6150". 
The  results  are  tabulated  In  Table  11.  In  view  of  Equation  (2),  the  response 
frequency  should  be  directly  proportional  to  the  square  of  the  diameter  if 
all  other  terms  are  held  constant. 

In  order  to  make  a  comparison  of  both  varying  length  and  varying  dia¬ 
meter,  theoretical  curves  are  plotted  with  extended  dimensional  decrements. 
These  are  shown  as  Figures  11  and  12.  On  Figure  12  a  linear  regression 
analysis  was  made  to  both  experimental  and  theoretical  curves.  The  slope 
of  the  experimental  curve  is  1169.6  and  the  slope  of  the  theoretical  curve 
is  1268.7.  This  agreement  appears  most  satisfactory  for  the  sample  size 
used.  Note  that  the  theoretical  cuirve  was  obtained  with  n  =  1.5186  same 
as  the  length  varying  case. 
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Figure  11. 

No-  Load  Round  Bars  Theoretical  Study 
(  J)=  0-7078  in.  L-  changes  ) 
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TABLE  II 


DIAMETER  CHAMGE  OF  ROUND  BAR 
LENGTH  12.037" 


Diameter 

In 


Weight  Frequency,  Hz 

Grams  -■  Exp  Theo 


7078 

605 

898 

389 

6995 

590 

888 

385 

6950 

585 

888 

381 

6880 

575 

894 

377 

6860 

565 

890 

378 

6800 

553 

887 

376 

6485 

500 

848 

360 

6250 

470 

816 

344 

6150 

455 

812 

339 

Rolled  Bars 


The  29  rolled  bars  received  from  Frankford  Arsenal  have  dimensions  of 
approximately  %  x  2  x  5^^  Inches.  Th&y  have  varying  hardness  ranging  from 
92.5  to  95  on  the  Rockwell  "B"  scale. 

The  experimental  set  up  was  consistent  In  that  a  sytrofoam  support 
apparatus  was  used  and  the  bars  were  struck  with  the  solenoid  used  In  pre¬ 
vious  experiments.  The  variability  of  the  resulting  frequencies  were  so 
small  that  no  difference  could  be  recognized  although  there  was  a  slight 
varying  In  hardness.  One  bar,  however,  gave  a  frequency  of  2  resolution 
width  below  the  others.  It  was  x-rayed  with  an  average  bar  of  the  medium 
hardness  of  the  range  group  but  no  flaws  was  detected  In  either  bar.  Fur¬ 
ther  examination  showed  that  the  bar  that  gave  the  lower  frequency  was  longer 
than  any  other  bars  by  0.1406  Inches.  According  to  the  theoretical  equation, 
length  Is  Inversely  proportional  to  frequency,  which  would  tend  to  explain 
the  lower  frequency. 


Thin  Wall,  Metal  Tubing 

Two  0.230  Inch  OD,  stainless  steel  tubes,  four  feet  long  were  received 
to  be  tested.  They  were  sent  by  WADCO  (a  subsidiary  of  Westlnghouse  Electric 
Corporation)  In  Richland,  Washington. 
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The  purpose  of  the  test  was  to  obtain  an  audio  range  envelope  by 
striking  the  tubes  at  different  points  along  the  tube  to  ascertain  any 
difference  In  frequency  that  could  be  attributable  to  some  defect  In 
the  tube. 

Many  different  conflguratlonr  were  tried,  two  of  which  are  shown  In 
Figures  13  and  14.  In  Figure  13,  the  tube  was  simply  supported  at  Its 
ends  by  two  pencils  as  shown.  It  was  then  excited  at  the  three  positions 
A,  B,  and  C.  The  envelope  was  recorded  and  played  back  to  the  hybrid 
computer.  The  results  are  shown  below  the  figure. 

In  Figure  14,  a  second  position  was  set  up  with  the  tube  being  sup¬ 
ported  in  a  sound  box  Is  shown.  Response  was  recorded  by  striking  the  tube 
at  points  a,  b,  c  and  d.  The  sound  was  then  played  back  to  the  computer 
with  the  results  as  shown  In  the  figure. 

Obtaining  an  envelope  from  the  tube  was  difficult  because  of  the  small 
mass  and  the  Inherent  flexibility  of  the  tube.  However,  the  vertical  . 
configuration  (Figure  14)  shows  promise  because  of  the  amplifying  power  of 
the  sound  box.  If  further  Investigation  could  be  made,  a  more  precise 
analysis  of  the  tubes  might  be  obtained. 

In  Summary,  further  Investigations  appear  worth  pursuing  but  since 
the  data  collecting  configuration  Is  constructed  for  machine  parts  of  some 
mass.  It  will  have  to  be  altered  to  obtain  satisfactory  response  from 
light  mass  tubing. 


Sensitivity  to  a  Small  Hole 


To  determine  the  sensitivity  of  the  acoustical  technique  to  discontinui¬ 
ties  a  1/60  Inch  hole  was  drilled  thru  the  center  diameter  axis  of  a  round 
carbon  steel  bar.  The  bar  was  tested  before  and  after  drilling  In  the 
usual  manner. 

The  first  run  to  ascertain  the  effect  of  the  1/60"  drilled  hole  through 
the  center  axis  of  the  round  carbon  steel  bar  was  not  successful.  The  bar 
had  a  diameter  of  0.70  Inch.  The  frequency  of  the  bar  before  the  hole  was 
drilled  was  within  the  Interval  of  982  -  986  Hz  which  Is  within  the  minimum 
frequency  resolution.  The  frequency  after  the  hole  was  drilled  remained 
within  the  982  -  986  Hz  Interval.  One  hypothesized  reason  was  that  the  un- 
drllled  bar  had  a  frequency  of  986  Hz  and  after  drilling  the  frequency  was 
only  as  low  as  982  Hz.  This  Is  one  of  the  fallacies  of  the  technique,  but 
the  probability  of  this  occurring  regularly  Is  not  great. 

The  second  run  was  more  successful.  The  diameter  of  the  bar  was  0.50 
Inch.  The  frequency  before  drilling  the  1/60"  hole  was  approximately  656  Hz. 
After  the  hole  was  drilled  the  frequency  dropped  to  approximately  640  cps. 
Note  the  difference  In  diameter. 
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since  the  theoretical  model  agrees  extremely  well  with  the  experi¬ 
mental  findings,  the  system  can  now  predict  the  frequency  change  due  to 
the  dimensional  changes.  As  a  result,  the  dimensional  changes  In  the 
stress  experiment  can  now  be  separated  from  the  overall  frequency  changes. 
That  Is,  the  frequency  change  due  to  stress  alone  can  be  Isolated. 

The  frequency  change  due  to  the  dlmens  *.onal  change  can  be  derived 
from  the  theoretical  equations.  The  relations  between  the  frequency  and 
the  length  or  diameter  are 

f  ~  1/L^  or  f  .s, 


Taking  a  derivative  of  the  above  equation  and  dividing  the  derivations 
by  the  original  equations: 


f 


If  there  Is  a  one  percent  change  In  length  (dL/L  =  1%)  or  In  diameter 
(dD/D  ■  IX)  a  two  percent  change  In  frequency  will  result.  Since  the 
frequency  change  due  to  the  dimensional  change  Is  Isolated,  the  stress 
and  frequency  relation  Is  established.  (Tile  minus  sign  means  an  Increase 
In  length  results  a  decrease  In  frequency). 


Conclusions --Phase  III 


Perhaps  the  most  exciting  and  unexpected  results  of  this  research 
was  the  high  degree  of  sensitivity  of  the  system  to  detect  differences 
In  stresses  In  machine  elements.  From  this  It  might  be  concluded  that 
the  possibility  might  exists^  for  determining  the  degree  of  stress  In 
machine  elements.  In  situ.  Elements,  such  as  pins,  connecting  rods, 
shafts,  and  linkages,  could  be  tested  without  being  removed  from  the 
system  and  any  unusual  degree  of  stress  In  those  members  could  be  as¬ 
certained.  Applications  of  this  type  of  nondestructive  test  could 
apply  to  such  diverse  systems  as  Internal  combustion  engines,  trans¬ 
mission  towers,  bridges,  and  airframes. 

As  V9B  anticipated,  the  system  had  the  ability  to  determine  changes 
In  dimensional  lengths  and  diameters  of  the  crder  of  magnitude  suffi¬ 
ciently  small  to  consider  the  process  here  developed  as  an  Inspection 
tool. 


For  many  years  the  control  of  liquids  and  gases  has  been  a  process 
which  could  easily  be  automated.  The  automated  control  of  hardware  Is 
In  the  Infant  stage  In  that  the  development  of  alrgages  and  other  auto¬ 
mated  Inspection  tools  are  new.  The  system  developed  here  could  add 
to  the  small  number  of  devices  now  available  to  automatically  Inspect 
machine  elements  or  any  other  matalllc  element  where  some  dimensional 
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property  is  of  Interest.  In  addition,  where  it  is  of  interest  to  deter¬ 
mine  significant  deterioration,  it  is  possible  to  apply  the  nondestructive 
test  to  determine  whether  the  element  is  in  the  essentially  original 
configuration. 

The  rolled  sections  received  from  Frankford  Arsenal  presented  a  rather 
unusual  situation.  These  items  of  steel  were  rolled  in  such  a  manner 
that  the  configuration  was  not  necessarily  uniform.  The  system  had 
the  ability  to  choose  one  of  29  bars  which  was  significantly  different 
in  configuration  from  the  others.  Since  the  bar  was  significantly  different 
(on  the  order  of  magnitude  of  0.25  inches),  the  sensitivity  of  the  system 
to  rolled  elements  remains  to  be  ascertained.  If  the  system  proves  applic¬ 
able  to  this  type  of  test,  then  the  whole  field  of  consideration  of  unusual 
size  and  shape  castings  might  be  amenable  to  this  type  of  nondestructive 
test. 


It  is  considered  that  the  sample  size  used  for  the  long,  thin  walled 
tubing  was  too  small  to  draw  any  meaningful  conclusions.  It  is  recognized 
that  there  are  mounting  and  damping  problems  which  tend  to  make  the  energy 
envelope,  received  from  vibrating  this  type  of  configuration,  difficult 
to  retrieve.  Nevertheless,  there  appears  to  be  no  real  barrier  to  devel¬ 
oping  a  computerized  nondestructive  testing  technique  for  these  types 
of  tubes. 

Again,  as  anticipated,  the  drilling  of  a  1/60  inch  round  hole  in  a 
machine  bar  gave  a  sufficiently  different  response  frequency  so  tha*-  any 
anticipated  nondestructive  computerized  test  to  pickup  faults  as  big 
as  1/60  inch  is  feasible.  Additional  work  will  be  necessary  to  refine 
the  tests  for  smaller  faults. 


Recommendations — Entire  Study 

Since  nearly  all  of  the  tests  considered  here  appear  to  be  feasible 
it  is  recommended  that  additional  work  within  each  area,  i.e.,  stress, 
dimensional  differences,  and  structural  integrity  be  investigated  fur¬ 
ther.  Standards  for  mounting  specimens  must  be  developed  and  background 
frequency  responses  for  standard  pieces  must  be  ascertained. 

The  sensitivity  of  all  of  the  tests  remains  to  be  investigated. 

It  Is  believed  by  the  researchers  that  a  small  portable  computer  can 
be  developed  to  conduct  the  tests  described  here.  ThiJ  would  require  the 
joint  effort  of  those  Interested  in  nondestructive  testing  and  those 
knowledgable  about  computer  and  circuit  design. 

In  summary,  the  overall  conclusion  from  this  research  indicates  that 
the  computer  will  assume  an  increasing  role  in  nondestructive  testing 
and  that  the  computer  docs.  Indeed,  have  the  ability  to  store  and  compare 
data  so  that  It  can  be  a  deci.>Lon  making  tool. 
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APPENDIX  A 


ANALYSIS  THEORY  AND  DATA  REDUCTION 

The  data  gathered  from  striking  a  bar  are  primarily  of  a 
transient  type  with  superimposed  resonance  modes  and  some  ran¬ 
dom  noises.  The  factors  that  contribute  to  the  shape  of  the 
energy  envelope  might  very  well  be  evolved  with  the  way  the 
bar  is  struck,  the  bar  size,  the  bar  material  and  the  manner  in 
which  the  system  is  mounted.  Two  approaches  seem  to  offer  the 
most  promise  in  analyzing  this  type  of  data.  The  direct  ap¬ 
proach  is  the  comparison  of  the  raw  data  of  the  defective 
sample  and  a  standard.  Another  approach  is  die  comparison  of 
the  data  in  their  frequency  domain,  which  is  obtained  by  a  trans¬ 
form  of  the  raw  data.  The  second  approach  was  chosen  for  the 
following  reasons:  1)  The  comparison  in  the  frequency  domain 
is  a  shorter  one,  especially  when  data  records  are  long.  2) 

The  amplitudes  in  the  frequency  domain  are  sensitive.  It  is 
anticipated  that  resonance  peaks  could  provide  easily  identifi¬ 
able  features  in  order  to  distinguish  the  defective  samples 
when  their  frequency  or  amplitudes  deviated  from  that  of  the 
standard  sample.  It  is  expected  that  a  resonant  frequency  will 
be  associated  with  each  sample  and  that  it  can  be  readily  iden¬ 
tified,  3)  Finally,  the  frequency  analysis  is  a  well  known 
and  tested  method  and  the  statistical  tools  are  readily  avail¬ 
able. 

The  transformation  from  the  time-history-records,  x(t),  to 
their  frequency  representation,  F(f),  is  generally  done  by  the 
Fourier  transform 


F(f) 


=  P(0 


exp(-2nift)  dt. 


Where  1  is  VT.  A  faster  version  of  the  transform  has  been 
developed  in  recent  years.  It  was  known  as  the  Cooley  and  Tukey 
Fast  Fourier  Transform.  The  analysis  program  developed  for  this 
research  project  has  adapted  a  Fast  Fourier  Transform  Subrou¬ 
tine  through  IBM  contributed  program  library  (subroutine  FOURl) . 
This  subroutine  performs  a  discrete  Fourier  transform  of  a  set 
of  N  equally  spaced,  time  samples,  x  ,  x.,  . . . ,  x^^  j  and  is 
given  by  ° 


N-1 

F  =  2  X  exp(-2TTij/N),  k  =  0,1,..., N. 
J=0  J 
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To  compute  the  frequency  spectrum  P,  ,  would  normally  require 
complex  multiplications.  However,  a  reduction  in  operations 
can  be  achieved  if  N  is  a  power  of  2.  An  algorithm  developed 
by  Cooley  and  Tukey^,  called  the  Fast  Fourier  Transform,  per¬ 
forms  the  same  transform  in  N(log2N)  complex  multiplications, 
thus  reduces  the  computing  time  by  approximately  . 

Mathematically  speaking,  the  trans  ormed  values,  F,  ,  are 
the  coefficients  of  the  Fourier  series  which  may  be  eimer  ne¬ 
gative,  posicive,  or  zero.  It  was  found  more  expedient  to  com¬ 
pare  the  power  spect-^al  density  function, that  is,  the  mean 
square  value,  rather  than  the  frequency  spectrum.  Apparently 
the  time  series  data  was  not  probabilistic  in  nature,  but  using 
the  average  squared  values  has  the  advantage  of  emphasizing  the 
peaks  and  de-emphasizing  the  minor  features.  For  this  reason 
the  power  spectral  density  function,  G  (f),  were  calculated 
which  were  obtained  by  filtering  the  sample-time-history  record 
in  a  frequency  range  from  f  to  (f  +  A  f) ,  squaring  them  and 
averaging  the  squared  output  from  the  filter.  Mathematically, 
the  power  spectral  density  function  is  defined  as^ 


G:.(f) 


lim  lim 
Af-*0  T— »  00 


1 

(Af)T 


,f,Af)  dt. 


where  T  is  the  total  record  time.  For  stationary  datA  ,  the 
power  spectral  density  function  is  related  to  the  auto-correla¬ 
tion  function  by  the  Fourier  transform 


G  (f) 
x' 


2  J  \(t)  exp(-2TiifT)  dr, 


where  R  (t)  is  the  auto-correlation  function  which  describes  the 
dependence  of  thesample-time-history  records,  x(t)  and  x(t  +  t), 
and  T  is  the  time  lag  between  the  two  records.  It  is  defined 
mathematically  as  an  average  over  the  observation  time,  T,  and 
is  the  real -valued  function. 


-  “I.  I 

Note  that  when  t  =  0,  the  auto-correlation  function  is  the  mean 
square  value  and  has  reached  a  maximum.  In  physically  realiz¬ 
able  situations,  only  one-sided  power  density  functions  can  be 
obtained,  where  frequency,  f,  varies  from  0  to  ».  That  means 
R  (t)  should  be  treated  as  an  even  function  of  t,  and  the  spec¬ 
tral  density  function  is  equivalent  to  folding  over  the  negative 
terms  from  the  Fourier  transform  into  the  positive  direction. 
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G  (f)  =  4  R  (t)  coe  Znfrdr 
X  Jq 


A  convenient  way  to  make  a  comparison  between  bars  Is  that  of  using 
the  cross  spectral  functions.  A  computer  program  has  been  written  to 
cross-correlate  the  data  between  bars.  Test  cases  have  been  tried  on  the 
program.  The  cross-correlation  function  is  a  function  which  describes  the 
mathematical  dependence  of  the  time  history  records  similar  to  the  auto¬ 
correlation  function  except  that  the  correlations  are  between  two  different 
sets  of  data.  It  can  be  expressed  as 


R 

xy 


(T) 


lim  1 
T  ^  oa  T 


1  x(t)y(t  +  T)  dt, 

"o 


which  is  a  real-valued  function  but  does  not  necessarily  have  a  maximum 
when  T  =  0,  nor  is  it  an  even  function.  The  cross-spectral  density  fun¬ 
ction,  Gjjy(f),  ie  obtained  by  the  Fourier  transform  of  the  cross -correlation 
function  Rjjy(T) .  It  is  done  in  a  smlliar  manner  as  the  auto-correlation 
function  Is^ transformed  into  the  power  spectral  density  function.  However, 
the  cross-correlation  function  is  not  an  even  function,  the  resulting  cross¬ 
spectrum  of  the  Fourier  Transform  has  the  complex  form 


G  (f)  =  C  (f)  -  iQ  (f), 
xy  xy  xy 

where  C^^yCf)  is  the  co-spectral  density  function,  which  is  defined  as 

'  4f““o  ^  dt. 

0 


It  represents  the  average  product  of  x(t)  and  y(t)  within  a  narrow  frequency 
interval  between  f  and  (f  +  At),  divided  by  the  frequency  interval.  The 
frequency  Qjjy(f),  is  called  the  quandrature  spectral  density  function.  It 
is  the  imaginary  part  of  the  transformed  result  and  is  defined  as 


0  i““.  J  ■it. 

0 
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where  y*(t,f,Af)  represents  the  time  shifted  y(t),  within  t'hfe  frequency 
Internal  f  and  (f  +  Af),  which  results  in  90  degree  phase  shift  in  the  com> 
plex  plane.  The  actual  computations  of  C^y  and  '  are  obtained  from  the 
transform  of  the  even  and  odd  parts  of  the  cross-correlation  functions,  de¬ 
noted  by  A  and  B  ,  respectively.  ,  '  i 

xy  xy  i  I 


A  (f)  =  [R  (f)  +  R  (f)1 
xy  2  xy  yx  -* 


I  ’  ' 

B  (f)  =  i  [R  (f)  -  R  (f)] 

I  xy  2  xy  yx  •*  '  .  ,  i 

I 

The  co-spectral  density  function  is  the  Fourier  transform  ofiA^^y  and  the 
quandratlve  spectral  density  function  is  the  Fourier  transform  of 
magnitude  of  the  cross  spectral  density  function  is  ^  •' 


'  G  (f)  =  (f)  +  (f),  , 

■  xy'  ,  xy'  '  xy'  *  > 

and  the  phase  angle  is 

‘  I 


I 

The  comparison  between  the  ^o  history  records  can  be  assessed  by  the 
cross  and  auto-spectral  density  functlpns.  The  combined  function  is  called 
the  coherent  function  and  is 'defined  as 


(f) 

xy^^^  “  G^(f)Gy(f)  ^ 


I 


;  2 

which  varies  between  0  and  1.  “  0»  indicates  incoherent  and  the 

value,  (f)  ^  indicates  fully  coherent.  If  two  bars  have  a  perfect 
match  thi^coherent  function  should  be  1  at  al,l  frequencies.  When  there  are 
deviations  between  bars,  the  coherence  function  drops  off  a  frequencies  where 
pdwer  spectra  are  out  of  agreement. 
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APPENDIX  B 

I  I  PROGRAMMING  AND  THEORY 

During  the  Initial  part  of  the  program  use  was  made  of  a 
digital  computer  (IBM  360,  Model  6!>)  entirely.  Then  analysis 
was  done  on  a  hybrid  computer  of  sufficient  capacity  to  analyze 
any  data  which  may  be  encountered  In  this  project.  The  hybrid 
system  In  question  uses  a  Sigma  5  for  the  digital  components 
and  EAI'  680  £or  the  analog  component. 

Since  hybrid  computers  are  relatively  new  Items  of  equipment 
It  may  be  well  here  to  mention  some  advantages  and  disadvantages 
of  analog,  digital  and  hybrid  computers. 

I 

The  analog  computation  has  the  disadvantage  of  having  a 
limited  precision  of  one  part  In  one  thousand.  Another  disadvan¬ 
tage  lies  In  the  necessity  to  Increase  the  size  of  the  computer 
capacity  directly  with  the  problem  size.  Further  It  Is  necessary 
to  establish  a  program  by  inner  connecting  computer  elements  with 
wires  and  the  reliability  limitations  Inherent  In  linear  devices 
(as  compared  to  discrete  devices). 

The  digital  computer,  on  the  other  hand.  Is  a  very  high  pre¬ 
cision  piece  of  equipment  with  nearly  absolute  repeatability. 
Further,  digital  programs  can  be  escaladed  In  complication,  a 
piece  at  a  time,  and  changing  from  one  program  to  another  Is 
simple,  iynbther  attractive  feature  of  the  digital  computer  lies 
In  Its  ability  to  handle  logical.  Boolean  and  algebraic  operations 
and  to  react  to  Inputs  that  send  It  on  completely  different  courses 
of  action,  such  as  Interrupts. 

Of  course,  the  digital  computer  has  Its  disadvantages,  the 
first  being  It  Is  a  slow  Integrating  device.  In  addition,  man- 
machine  Interface  continues  to  be  a  time  consuming  problem. 

Finally,  cost  are  very  great  In  digital  equipment. 

The  main  advantage  of  the  analog  computer  lies  In  Its  ability 
to  receive  voltage  equivalent  signals  and  deal  with  these  directly, 
whereas  the  Information  must  be  digitized  to  employ  a  digital 
computer . 

The  hybrid  computer  simply  combines  the  Integrating  speed  of 
the  analog  and  the  acci  racy,  arithmetic  and  logical  operation  of 
th3  digital  so  that  problems  done  either  poorly  or  not  at  all  by 
either  analog  or  digital  machines  can  be  solved. 
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For  the  particular  application  of  hybrid  computer  to  the  in¬ 
vestigation  in  question  the  analog  to  digital  converter  and  the 
digital  part  of  the  computer  allow  the  on-line,  real-time  anal¬ 
ysis  of  the  data  so  that  information  can  be  obtained  immediately. 

The  sound  signature  analysis  is  usually  done  by  converting 
the  time  series  data  into  their  frequency  domain  by  a  direct 
Fourier  intergal  transformation  of  the  data.  The  resulting  fre¬ 
quency  spectrum  describes  the  frequency  response  of  the  sounding 
system.  The  frequency  representation  of  a  sound  signature  re¬ 
veals  much  of  the  physical  characteristics  of  the  sounding  source 
When  the  relative  power  representation  of  the  sound  signature  is 
described  as  a  function  of  frequency,  a  power  spectral  density 
function  is  obtained  indirectly  by  the  Fourier  transform  of  its 
autocorrelation  function.  Hie  distinction  between  the  power 
spectrum  and  the  direct  Fourier  transformed  spectrum  from  the 
sound  signature  is  similar  to  the  comparison  of  a  sine-squared 
function  to  a  sine  function.  The  natural  frequency  is  greatly 
magnified  and  the  noise  frequencies  are  relatively  depressed  in 
the  power  spectral  representation.  No  negative  values  should  be 
found  in  the  power  spectrum  theoretically  (although  digitized 
formulation  sometimes  produces  a  slight  negative  spike).  The 
reason  that  the  power  spectral  analysis  is  favored  lies  in  the 
ascertaining  of  whether  a  minute  change  in  the  specimen  can  cause 
corresponding  change  in  the  magnified  response  frequency  peak  of 
the  power  spectrum. 

To  illustrate  the  direct  Fourier  transform  of  a  time  series 
data  x(t),  let  F(f)  represent  their  transformed  value  in  the  fre¬ 
quency  domain. 


F(f)  -  r  x(t)’e‘^^^^*^dt 

where  f  is  the  frequency  variable,  i  is  equal  to  \J-1,  and  t  is 
the  time  variable.  If  x(t)  is  a  single  rectangular  pulse,  see 
Figure  15,  the  corresponding  frequency  spectrum  is  shown  in  Figure 
16. 


Its  power  spectrum  representative  is  shown  in  Figure  17. 

If  x(t)  is  a  pure  sine  function,  its  corresponding  frequency 
representation  is  shown  in  Figures  18  &  19.  The  sound  signature 
in  the  heat-treated  bar  experiments  were  sinusoidal  waves  superim 
posed  with  an  exponential  damping,  with  some  low  frequency  com¬ 
ponents  and  high  frequency  noise.  The  recorded  sound  data  were 
first  digitized  by  either  a  CAT  computer  or  hybrid  system.  Its 
autocorrelation  functio:i  was  determined  by  che  relation. 
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Frequency 


x(t)  x(t  +  t)  dt, 


R(t) 


where  t  represents  the  time  lag  between  the  sampling  points  and 
T  Is  the  total  sample  record  length  in  time.  The  time  lag,  t,  is 
a  variable  ranging  from  0  to  the  maximum  lag,  L.  To  get  a  small 
statistical  uncertainty  choose  small  L,  such  that  L  <<  T.  But 
in  order  to  have  high  frequency  resolution  choose  a  larger  L, 
such  that  L  <  T.  A  compromise  choice  was  used  in  the  sound  data 
analysis  such  that  L/T  »  ^  to  %  where  T  represents  1024  sampling 
points. 


The  time  Interval  between  sampling  points  could  be  selected 
using  the  CAT  unit.  The  most  frequently  used  time  interval  was 
0.25  milllsec.  This  would  give  a  Nyqulst  frequency  f^  (cut  off 
frequency)  of 


‘c  -  i  “  S  ■  20“ 

where  h  is  the  time  Interval  in  seconds. 

For  the  hybrid  system,  the  time  interval  between  points  is 
controlled  by  a  digitizing  circuit.  The  perferred  interval  is 
0.2  mlllisec,  which  gives  a  Nyquist  frequency  of  2500  Hertz. 

After  the  digitized  information  was  obtained,  it  was  trans* 
ferred  and  stored  in  the  memory  core  of  the  digital  computer.  The 
autocorrelation  function  R(jl)  is  then  calculated  by  the  digitized 
formula 

N-I 

'v«.^ 

n-i 

where  Jl  is  the  lag  number,  L  is  the  maximum  lag  number  and  R(.^) 
is  the  autocorrelation  function  at  lag  X  •  '^he  equivalent  band¬ 
width  of  the  computational  procedure  is 


where  the  equivalent  bandwidth  B  is  twice  the  range  found  by  di¬ 
viding  the  frequency  Interval  (0,  f  )  into  L  equally  spaced  parts, 

or  f  /L  parts.  ® 

c 


0,1,2, .  .  .L 
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The  auiL^/Correlatlr-n  function  may  take  on  negative  as  well 
as  positive  values.  When  i  “  0,  the  autocorrelation  function 
is 

N 

■  “'<»  -  F  I  <”n>"  - 
n=l 

where  x  are  the  digitized  zero-mean  time  series  data.  R  is  re¬ 
lated  to  the  sample  variance  s^  by  the  relation 


R 

o 


/N^IN 

V  N ; 


2 

s 


The  power  spectral  density  function,  G(f),  is  calculated  by  means 
of  a  fast  Fourier  transform  algorithm  developed  by  Cooley  and 
Tukey,  which  reduces  a  complex  multiplications  to  only  L  lag^L 
complex  multiples  it  the  max  lag  L  is  an  Integer  power  of  2, 

This  reduces  the  computation  time  by  as  much  as  99%.  The 
definition  of  power  spectral  density  function  is  given 
in  Appendix  A.  However,  for  continuity  of  dis¬ 
cussion,  some  repetition  will  be  done  here.  G(f)  is  always  a 
roal”vaiued,  non-negative  function.  By  definition 


G(£)  «  lira  lira  J  dt 

^£-*“0  T*^®  o 

where  x(t,£,A£)  is  that  portion  o£  x(t)  in  the  frequency  range 
from  f  to  (£fA£)«  For  stationary  data,  the  power  spectral  density 
function  and  the  autocorxelation  function  is  related  by  a  Fourier 
transform  as  follows: 

G(f)  -  2J  R(T)e"^’^^^  dT 

Since  R(t)  Is  sn  even  function  of  t,  its  digital  relation  reduces 
to 

L-1 

. 

X-0 

where  W  ■  e  L  ts  the  max  lag,  Rj^  is  the  digitized  auto¬ 

correlation  function  ar.i  G^^  is  the  digitized  power  spectrum. 
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Appendix  c 


Stressed  Bar  Model  Theory 


Applying  d'Alasobert's  principle  to  a  fre'i  lateral  vibration  of  ( 
supported  bar  we  have: 


El  24  -  s  si 


ax 


_  .  0  ^ 
2  p  2 
ax  at 


where  E  =  Young's  modulus  (psl)  '  1 

4 

I  =  moment  of  Inertia  of  crofs  section  area  (In  ) 
p  =  linear  density  (slug/ln) 

S  =  applied  force  (#) 

I  ' 

t  =  time  (sec) 

I  ..II 

Boundary  conditions:  y(L,t)  <>  0 

,  I  y(0,t)  -  0  ' 

I 

Propdsed  solution:  ' 


-I 


A  (t)  Sin  + 
n  .  L 


n  «  0 


I 


B„(t)  Co,  ms 


n  •=  0 


B^(t)  *=  0  due  to  the  boundary  condition 


ILI 


ax 


n  »  0 


L 


I  I 


2  »  2  2 

M  =  -  y  ^  Sin  ^ 

ax^  n  -  0  ^ 


f  i 


52 


simply 


/ 


I 


I 


I 


H  ■  I  ^ 

n  =  0 


Substituting  back  into  Equation  (2) 


4  4 


^  El  ill  AJt)  Sin  ^  S  A__(t)  Sin  !f- 


2„2 

n  n 


n  =  0 


n  «  0 


■I 

n  =  ,0 


P  A„(t)  Sin  ^ 


The  eigen  Equation  is  then: 


4  4  2  2 

El  A^(t)  +  S  A^(t)  +  p  A^(t) 


or,  simply  a  wave  equation 


^n(«  ^  [  i  ^  ^  ^  ■  “ 


The  solution  of  the  wave  equation  is 


A  (t)  «=  C  Sin  cot  +  r  Cos  cut 
n  n  n 


with  arbitrary  coefficients  C  and  D  . 

n  n 
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The  correspoodirig  frequency  li 


or,  the  response  frequency  is: 
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APPENDIX  D 


1040  itressed  Bar  Data 


Load 

lb 

Depth 

In 

Width 

in 

Length 

in 

Frequency, 

Exp 

Hz 

Thfco 

.2534 

.9994 

11.990 

242 

175 

2000 

.2534 

.9994 

11.993 

367 

197 

3000 

.2534 

.9994 

11.999 

382 

216 

Q> 

4000 

.2534 

.9994 

12.002 

398 

234 

n 

5000 

.2533 

.9994 

12.005 

414 

251 

.2533 

.9994 

12.007 

429 

266 

7000 

.2533 

.9994 

12.010 

— 

281 

8000 

.2533 

.9994 

12.012 

453 

295 

9000 

.2532 

.9994 

12.015 

— 

309 

10000 

.2532 

.9994 

12.017 

— 

321 

1000 

.2528 

12.340 

351 

164 

2000 

.2528 

12.342 

351 

191 

3000 

.2528 

12.344 

375 

210 

w 

4000 

.2528 

12.347 

332 

228 

Q> 

5000 

.2528 

12.352 

398 

245 

6000 

.2528 

12.357 

414 

260 

to 

6900 

.2527 

12.362 

988 

273 

8000 

.2527 

12.365 

960 

289 

8900 

.2526 

12.368 

945 

301 

10000 

.2526 

12.371 

930 

315 

Inches 


Bar 

A 

B 

1 

.811 

.965 

2 

.692 

.718 
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Labels  A  and  B  indicate  the  length  of  the  bar  was  used  for  the 
gripping.  These  values  will  be  given  below  each  remaining  group 
of  tables. 
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Load 

Depth 

Width 

Length 

Frequency, 

Hz 

lb 

ill 

In 

in 

Exp 

Theo 

1000 

.2534 

1.000 

12.279 

335 

170 

2100 

.2534 

1.000 

12.285 

281 

194 

3000 

.2534 

1.000 

12.290 

250 

212 

R) 

4000 

.2534 

1.000 

12.294 

226 

230 

0) 

5000 

.2533 

1.000 

12.298 

195 

246 

6000 

.2533 

1.000 

12.302 

414 

261 

w 

7000 

.2533 

1.000 

12.305 

429 

276 

8000 

.2533 

1.000 

12.309 

953 

290 

.2532 

1.000 

12.314 

929 

303 

10000 

.2532 

1.000 

12.317 

— 

316 

1000 

.2540 

.9975 

11.952 

351 

177 

2000 

.2540 

.9975 

11.955 

375 

199 

3000 

.2540 

.9975 

11.958 

390 

219 

a 

4000 

.2539 

.9975 

11,962 

406 

236 

5000 

.2539 

.9975 

11.965 

414 

253 

6000 

.2538 

.9975 

11.969 

269 

7000 

.2538 

.9975 

11.972 

... 

284 

7900 

.2538 

.9975 

11.975 

... 

296 

9000 

.2538 

.9975 

11.980 

— 

312 

10000 

.2538 

.9975 

11.984 

••• 

324 

Inches 

Bar 

A 

B 

3 

.881 

.864 

4 

.945 

1.132 
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Load 

Depth 

Width 

Length 

Frequency, 

Hz 

lb 

in 

In 

in 

Exp 

Theo 

1000 

.2527 

1.020 

11.978 

343 

177 

2000 

.2527 

1.020 

11.983 

375 

199 

3000 

.2527 

1.020 

11.988 

'  390 

219 

3900 

.2526 

1.020 

11.993 

398 

235 

5100 

.2526 

1.020 

11.996 

414 

255 

*'v»8 

6000 

.2526 

1.020 

12.000 

429 

269 

7000 

.2526 

1.020 

12.005 

445 

283 

8000 

.2526 

1.020 

12.008 

921 

297 

9000 

.2525 

1.020 

12.011 

— 

311 

10000 

.2525 

1.020 

12.014 

--- 

324 

_ _ _ _ _ _ _ _ _ _ _ 1 

1000 

.2527 

.9998 

IBl 

.2527 

.9998 

.2527 

.9998 

12.242 

375 

213 

u 

.2526 

.9998 

12.247 

390 

231 

.2526 

.9998 

12.252 

185 

248 

o\ 

6000 

.2526 

.9998 

12.255 

414 

264 

7000 

.2526 

.9998 

12.258 

429 

278 

8000 

.2526 

.9998 

12.261 

445 

292 

9000 

.2525 

.9998 

12.265 

453 

305 

10000 

.2525 

.9998 

12.271 

914 

318 

Inches 

Bar 

A 

B 

5 

.956 

1.108 

6 

.780 

1.350 
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Load 

lb 

Depth 

in 

Width 

in 

Length 

in 

Frequency, 

Exp 

Hz 

Theo 

1000 

.2535 

1.035 

12.180 

312 

174 

2000 

.2535 

1.035 

12.183 

265 

196 

3000 

.2535 

1.035 

12.189 

375 

215 

4000 

.2534 

1.035 

12.196 

390 

232 

5000 

.2534 

1.035 

12.200 

406 

249 

6000 

.2534 

1.035 

12.205 

421 

264 

7000 

.2534 

1.035 

12.210 

429 

279 

7900 

.2534 

1.035 

12.213 

445 

292 

9000 

.2534 

1.035 

12.217 

— 

306 

10000 

.2533 

1.035 

12.221 

— 

319 

.2531 

1.005 

12.226 

320 

171 

.2531 

1.005 

12.230 

265 

193 

.2531 

1.005 

12.234 

382 

^12 

4000 

.2530 

1.005 

12.238 

390 

230 

5000 

.2530 

1,005 

12.241 

206 

247 

6000 

.2530 

1.005 

12.245 

421 

262 

7000 

.2530 

1.005 

12.248 

429 

111 

8000 

.2529 

1.005 

12.252 

445 

291 

9100 

.2529 

1.005 

12.255 

305 

10000 

.2529 

1.005 

12.260 

906 

317 

Inches 


Bar 

A 

B 

7 

.860 

.970 

8 

.815 

1.020 
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Load 

Depth 

Width 

Length 

Frequency, 

Hz 

lb 

In 

In 

in 

Exp 

Theo 

1000 

.2535 

.9988 

12.130 

296 

174 

2000 

.2535 

.9988 

12.134 

367 

195 

3000 

.2535 

.9988 

12.139 

382 

215 

tjj 

4000 

.2534 

.9988 

12.142 

398 

233 

Q> 

5000 

.2534 

.9988 

12.147 

414 

250 

=#: 

6000 

.2534 

.9988 

12.152 

— 

265 

VO 

7000 

.2533 

.9988 

12.156 

437 

280 

8000 

.2533 

.9988 

12.158 

445 

294 

9000 

.2533 

.9988 

12.161 

— 

308 

10000 

.2533 

.9988 

12.166 

898 

320 

1000 

.2528 

.9800 

12.123 

281 

172 

2000 

.2528 

.9800 

12.125 

367 

194 

ta 

3000 

.2528 

.9800 

12.128 

382 

214 

CD 

►I 

4000 

.2527 

.9800 

12.131 

398 

232 

5000 

.2527 

.9800 

12.135 

406 

249 

h-* 

O 

6000 

.2527 

.9800 

12.139 

421 

264 

7000 

.2526 

.9800 

12.142 

437 

279 

8000 

.2526 

.9800 

12.145 

937 

293 

9000 

.2526 

.9800 

12.150 

914 

307 

10000 

.2526 

.9800 

12.155 

898 

320 

Inches 

Bar 

A 

B 

9 

.880 

1.150 

10 

.922 

.970. 
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mm 


Load 

Depth 

Width 

Length 

Frequency, 

Hz 

lb 

In 

in 

in 

Exp 

Theo 

1|!||H 

.2526 

.9945 

12.082 

281 

173 

s 

.2526 

.9945 

12.085 

367 

195 

S 

.2526 

.9945 

12.090 

382 

214 

BR  S 

.2525 

.9945 

12.095 

398 

232 

Di 

t 

.2525 

.9945 

12.100 

414 

249 

Im 

.2525 

.9945 

12.105 

421 

264 

.2525 

.9945 

12.109 

437 

279 

t 

.2524 

.9945 

12.113 

929 

293 

BB  [ 

.2524 

.9945 

12.116 

910 

306 

10000 

.2524 

.9945 

12.020 

890 

321 

.2535 

.9965 

12.125 

343 

173 

2000 

.2535 

.9965 

12.130 

367 

194 

3000 

.2535 

.9965 

12.135 

382 

214 

dcJ 

4000 

.2535 

.9965 

12.140 

398 

323 

5000 

.2534 

.9965 

12.144 

414 

248 

=#= 

6000 

.2534 

.9965 

12.148 

421 

264 

N> 

1 

.2534 

.9965 

12.152 

437 

279 

.2534 

.9965 

12.155 

929 

293 

.2534 

.9965 

12.159 

914 

306 

10000 

.2533 

.9965 

12.164 

890 

318 

Inches 

Bar 

A 

B 

11 

.955 

1.00 

12 

.905 

1.00 
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Load  Depth  Width  Length  Frequency,  Hz 

lb  in  in  in  Exp  Theo 


1000 

.2525 

1.005 

12.120 

304 

173 

2000 

.2525 

1.005 

12.126 

367 

194 

3000 

.2525 

1.005 

12.131 

382 

214 

u 

4000 

.2525 

1.005 

12.135 

398 

232 

to 

H 

5000 

.2524 

1.005 

12.140 

414 

248 

=Sf! 

6000 

.2524 

1.005 

12.143 

421 

264 

(jJ 

7000 

.2524 

1.005 

12.146 

437 

279 

8000 

.2524 

1.005 

12.150 

929 

293 

9000 

.2524 

1.005 

12.155 

914 

306 

10000 

.2524 

1.005 

12.159 

898 

319 

1700 

.2521 

.9958 

11.505 

375 

198 

2050 

.2521 

.9958 

11.507 

390 

206 

3000 

.2520 

.9958 

11.514 

405 

224 

w 

a 

4000 

.2520 

.9958 

11.520 

420 

242 

n 

4900 

.2520 

.9958 

11.523 

--- 

257 

6000 

.2520 

.9958 

11.528 

452 

275 

7000 

.2520 

.9958 

11.532 

468 

290 

8000 

.2520 

.9958 

11.536 

— 

304 

9000 

.2520 

.9958 

11.538 

684 

317 

10000 

.2520 

.9958 

11.540 

700 

330 

Bar 

A 

B 

13 

.818 

1.105 

14 

NA. 

NA. 
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APPENDIX  E 


Theoretical  Model  to  Develop  Frequency  of 
Non-stressed,  Simply  Supported  Bars 


Simply  Supported  Bar  Without  Stress 

The  moment  equation  of  a  simply  supported  bar  is  given  as 

El  =  -M 
dx 

If  we  take  the  second  derivative  of  the  above  equation  with  respect 
to  X,  we  obtain  the  shearing  force 

(El  ■^)  =  F 
dx  dx 

now,  applying  d'Alembert's  principle,  the  vibrating  bar  is  loaded  by 
inertia  forces  varying  along  the  length  of  the  bar; 

■  '’72 

dt 

and  is  balanced  by  the  shearing  force 

dx  at 

where  E  =  Young's  modulus  (psi) 

p  -  linear  density  (slug/ln) 
y  =  y(x,t) 

4 

1  =  moment  of  inertia  of  cross  section  area  (in  ) 

Boundary  conditions:  y(L,t)  «=  y(0,t)  ■  0 

00  00 
Propose  solution:  y  “  ^  ^ 


n>0 


n*0 


B^(t)  »  0  due  to  the  boundary  condition. 
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/ 


I 


Sin  Sin  ; 


n“0 


at 


I  Stn  ^ 


n=0 


A.  I  sin®  ; 


.^To  '■■' 


Substitute  back  into  eqtiatlbn  (1). 


'4  4 


nnrx 

L 


P  ' 


n«0 


I 

1^0  ' 


Sin 


Eqtiating  coefficient  of  sine  functions: 
4  4 

El  ^  A^(t)  +  p  A^Ct)  -  0 


or, 


n<«  *  (j  -  0 


This  is 'a  wave  equation  whose  frequency  is 


(U 


yPT  4  4  2  2 

El  nrt_  .  nrr 

P  IF  IF 


The  solution  of  the  waye  equation  is: 

( 

A  (t)  <■  C  Sin  (ut  +  D  Cos  cut 
n  n  n 

The  response  frequency  is: 

f 


cu  ^  n^n  ,  /  El 

^  ay  p 


I  i 


»  I 
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I 
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